model;
{
b.clu[1:na] ~ car.normal(mapa[],weia[],numa[],tau.clu);

b.tem[1:nt] ~ car.normal(mapt[],weit[],numt[],tau.tem)
for (j in 1:nt) {eff.tem[j] <- b.tem[j]-b.tem[1]}

for (i in 1:na) {
for (j in 1:nt) {
ind[i,j]        ~  dnorm(media[i,j],tau.indi);
media[i,j]  <- b.ete2[i]+eff.tem[j];
                }
b.ete2[i]        ~  dnorm(0.00,tau.ete2);
                }

for (i in 1:na) {
for (j in 2:2){
oss[i,j]        ~  dpois(lambda[i,j]);
log(lambda[i,j])  <- log(att[i,j])+alpha+b.clu[i]+b.ete[i]+beta*ind[i,j+3];
rr[i,j]           <- lambda[i,j]/att[i,j];

}
b.ete[i]        ~  dnorm(0.00,tau.ete);
}


beta ~ dnorm(0.00,1.0E-3);

alpha           ~ dflat();
tau.ete         ~ dgamma(1.0E-3,1.0E-3);
tau.indi         ~ dgamma(1.0E-3,1.0E-3); 
tau.clu         ~ dgamma(nu.clu,sf.clu);
tau.tem         ~ dgamma(nu.tem,sf.tem);
tau.ete2        ~ dgamma(1.0E-3,1.0E-3); 
}